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Abstract 

Using 3.08 million ip(2S) decays observed in e + e~ collisions by the CLEO detector, we present 
the results of a study of the ip(2S) decaying into baryon-antibaryon final states. We report the 
most precise measurements of the following eight modes: pp, AA, (first observation), 

£+£+ (first observation), and E°S°, and place upper limits for the modes, and 
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The study of ip{2S) production in e + e~ and its subsequent decay into two hadrons 
provides a test of the predictive power of QCD 1] , including information on gluon spin, quark 
distribution amplitudes in baryon-antibaryon pairs, and total hadron helicity conservation. 
The decays are predicted to proceed via the annihilation of the constituent cc into three 
gluons or a virtual photon. This model leads to the prediction that the ratio of the branching 
fraction into a specific final state to the branching fraction of the J/\I/ into that same 
state should be a constant value of approximately 13%, the corresponding ratio for the 
dilepton final state [2]. This rule, which was previouslvrefereed to as the "12% rule", is 
roughly obeyed for several channels, but fails for others [3J. This Letter concentrates on the 
investigation of the two-body decays of the ip(2S) into a baryon-antibaryon pair. Previous 
measurements |H of the ip(2S) decaying into these final states have very large statistical 
uncertainties. The CLEO |5j detector, with the advantage of good charged and neutral 
particle detection efficiencies, together with secondary and tertiary vertex reconstruction, 
gives us the opportunity to make the most precise measurements yet of these branching 
fractions. 

The data used in this analysis were collected at the CESR e + e~ storage ring, which 
has been reconfigured to run near cc threshold by inserting 6 wigglers Our analysis is 
based on 3.08 xlO 6 ip{2S) decays which corresponds to the total integrated luminosity of 
5.63 pb -1 . Approximately half of the ip(2S) data (2.74 pb -1 ) were taken with the CLEO III 
detector configuration p and the remainder (2.89 pb -1 ) of the ip(2S) data, and all the 
continuum data (20.70 pb -1 , yfs = 3.67 GeV) were taken with the reconfigured CLEO-c 
detector We generated 10,000 Monte Carlo events using simulations of each of the two 
detector configurations for each of the eight decay modes, using a GEANT-based 0] detector 
modeling program. For all modes we generated Monte Carlo samples with a flat distribution 
in cos0, where 9 is the angle of the decay products in the center-of-mass system (in colliding 
beam experiments this angle is measured relative to the beam axis). For spin 1/2 baryons 
in the baryon octet, we then weighted the Monte Carlo samples with a 1 + cos 2 8 angular 
distribution, in agreement with the naive expectation. Possible deviations from these angular 
distributions will be one source of systematic uncertainty in our measurements. 

We begin by reconstructing the hyperons in the following decay modes (branching frac- 
tions 4] are listed in parentheses): A->p7r _ (63.9%), S + ^j9vr°(51.6%), S°^A7(100.0%), 
5-^Att-(99.9%), H°^A7r°(99.5%), S* ^S^tt+(51.6%), and fi-->Atf-(67.8%). To dis- 
criminate between protons, kaons, pions, and electrons, we combined specific ionization 
(dE/dx) measured in the drift chamber [^| and log-likelihoods obtained from the RICH sub- 
detector to form a joint log-likelihood difference: C(p — ir) = Lrich(p) — Lrich{^) + 
a dE/dx(p) ~ VdE/dxi 71 )' where the more negative C(p — ir), the higher the likelihood that 
the particle is a proton compared to a pion. Our requirement on these quantities varies in 
value from mode to mode depending upon background considerations. Further details of 
this procedure may be found elsewhere |8j. For protons in the pp and S° decay modes we 
require C(p — ir) < — 9 and C{p — K) < — 9 (3 o separation), and for those in the pp mode 
we require additional criteria C(p — e) < —9 as, in this case alone, electrons are a potentially 
significant background. For protons that are the daughters of A decays we make the looser 
requirements of C(p — ix) < and C(p — K) < 0, which are both very efficient. Kaons from 
Cl~ decays are strongly identified with C(K — n) < — 9 and L{K — p) < —9. Pions are 
only required to have energy loss measurements consistent with the pion identity. Photon 
candidates were identified in the Csl calorimeter. 

The analysis procedure for reconstructing A and S~ closely follows that presented ear- 
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Invariant Mass (GeV/c 2 ) 

FIG. 1: Inclusive hyperon invariant mass distribution in tp{2S) data: (a) A^p7r~, (b) X + — >p7r°, 
(c) S°^ 7 A, (d) S-^Avr", (e) S^Avr , (f) H*°^S-vr+, (g) n-->AK~. 

lier ||. The ft~ reconstruction follows the steps of the S~ reconstruction by replacing a 
pion by a kaon. S° hyperons are reconstructed by combining a A with a photon and requir- 
ing the photon energy in the crystal calorimeter to be in excess of 50 MeV, not matched 
to a charged track, and consistent in shape to that from a photon. S° and S + hyperon 
reconstruction is complicated by the fact that there is no direction information from the Csl 
photon clusters for the 7r° reconstruction. A kinematic fit is made to the hypothesis that 
the parent hyperon started at the beamspot, and decayed after a positive pathlength to a 
decay vertex which is the origin of the 7r° — ► 77 decay. A cut was placed on the \ 2 °f the fit 
to this topology, which includes the fit to the ir° mass from the newly found hyperon decay 
vertex. The S*° candidate is reconstructed by vertexing a positively charged pion candi- 
date, consistent with coming from the beamspot, with the already found H~. All hyperon 
candidates were required to have vertices significantly separated from the beamspot, with 
the flight distance of the A measured from the parent (S°, S~, or Q~) decay vertex, to be 
positive. In Figure ^ we show the inclusive hyperon yields in the ip{2S) data. We observed 
48xl0 3 A candidate events in the ip(2S) data as shown in Fig^a) with a fitted width (a) 
found to be 1.53 MeV/c 2 , consistent with Monte Carlo estimates. The fitted yields (widths) 
in Fig. ^b) and (c) for S + and E° candidates are approximately 3xl0 3 (4.0 MeV/c 2 ) and 
1x10 s (4.7 MeV/c 2 ), respectively. Invariant mass distributions for the S~ and S° are shown 
in Fig. ^d) and (e), respectively. The fit yielded approximately 2.7xl0 3 and 1.2xl0 3 for 
the signal events and 2.4 MeV/c 2 and 3.6 MeV/c 2 for widths for the S~ and S° candidates, 
respectively. The fits in Fig. Dp) and (g) yielded 115(4.69 MeV/c 2 ) and 47(2.69 MeV/c 2 ) 
events (width) for the S*° and Q~ spin 3/2 baryons. Hyperon candidates within 3a of their 
nominal masses are considered for further analysis. We then combined the 4-momenta of 
these baryons with the corresponding 4-momenta of their charge conjugates and formed 
ip(2S) candidates. 

We note that baryon-antibaryon pairs produced by ip(2S) decays have a back-to-back 
topology, and each decaying baryon has exactly the beam energy. Momentum conservation 
is imposed on the reconstructed baryon-antibaryon pairs by demanding the vector sum of 
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the total momentum in an event |Sp|/i?beam be less than 0.04. This eliminates background, 
via ttttJ/iP or 7Xcn (n = 0,1, and 2), which have extra tracks or showers. 

For each baryon-antibaryon candidate event, we calculated the scaled visible energy, 
Evis/ y/s, where E vis is the energy observed in an event and y/s is the center of mass energy. 
We define our signal region be 0.98 < E vis /\/s < 1.02, and two sideband regions of 0.94-0.98 
and 1.02-1.06 as representative of the combinatorial background. 

We also studied the continuum data to check for a possible contribution to our ip(2S) 
signals from this source. This was found to be non-negligible for four of the decay modes as 
shown in Table UJ We multiplied the yield from the continuum data by a scaling factor which 
was calculated taking into account the differences in luminosity (0.2720), a 1/s 5 correction 
(0.9572) for baryons Q, and the values of the efficiencies in the CLEO III and CLEO-c 
detector configurations before subtracting it from the ip{2S) yields. This scaling factor was 
0.2547 for the pp case and similar for the other modes. 

Other possible background sources are the daughters from the J j ip decay combined with 
the charged or neutral transition pions which can produce possible cross-feed in the signal re- 
gion. To estimate this cross-feed we generated corresponding Monte Carlo samples of 20,000 
events and looked for events which passed our selection criteria for the modes concerned. 

Figure El shows the scaled energy distribution for each of the decay modes. In all cases, 
clear signals are seen, with widths as expected from Monte Carlo simulation, and very little 
combinatorial background in the sideband regions. In each mode we calculate the number 
of ip(2S) decays to each final state ("signal yield") N s = S^S) ~~ A/)(2S) — S*f — fs 4 B c 
(where S^s) is the total number of events in the signal region, B^s) events in the scaled 
sidebands, S x f is the contribution from the cross-feeds, and f s ■ B c is the scaled continuum 
contribution with scaled sidebands subtracted) . These yields are shown in Table UJ The effi- 
ciencies, calculated using a weighted average of results from CLEO III and CLEO-c detector 
simulations, are also shown. 

We evaluated the following systematic uncertainties to our measured branching fractions: 
3% uncertainty on the number of ip(2S) decays in our sample; 1% uncertainty in the simu- 
lation of our hardware trigger; 1% uncertainty in the reconstruction in each charged track 
in the event; 1% uncertainty for proton identification of tracks coming from the beamspot, 
and 2% for proton and kaon identification of tracks coming from the secondary vertices; 
1% and 2% uncertainties for photon detection and ir° reconstruction, respectively; and 1% 
for background subtraction in AA mode. The detection efficiency also depends upon the 
angular distribution of the hyperons. We find a 10% change in efficiency when we change 
from a flat distribution in cos# to one of the form 1 + Acos 2 # (where A = 1). The value 
of A may not be strictly 1 because of effects due to finite charm quark mass J9| and hadron 
mass effect from 0(v 2 ) and higher twist corrections to the effective QCD Lagrangian. We 
compare the detection efficiencies for spin one-half baryon generated samples as discussed 
above with that obtained by the E835 collaboration |l0j and assign the difference, 3.3%, as 
the systematic uncertainty. We assign a 10% uncertainty in angular distribution to spin 3/2 
baryons, which covers all reasonable possibilities. The systematic uncertainty in the hyperon 
efficiency was estimated from a comparison of the efficiency of the various requirements in a 
data and Monte Carlo sample. This is largest for the S + and S° modes where the agreement 
is less satisfactory. In all modes the uncertainty in the reconstruction is doubled as there 
are two hyperons per event. 

Table |H] displays the breakdown of our systematic uncertainties. The last column tabu- 
lates the total systematic uncertainty added in quadrature for each mode. 
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FIG. 2: Scaled energy (E v i s /y/s) distribution for eight i^>{2S) decay modes: (a) pp, (b) AA, (c) 
£+£+, (d) (e) E~E~, (f) (g) E°ZP : and (h) n~TF. Signal Monte Carlo, data, and 

scaled continuum plots are shown in solid line, error bar, and filled histograms, respectively. 



TABLE I: Branching ratios of ip{2S) decaying to baryon-antibaryon pairs. The second, third, 
fourth, and fifth columns show the yields in the signal region, scaled sidebands, scaled continuum, 
and cross-feeds, respectively. The sixth and seventh columns tabulate the efficiencies and the 
measured branching fractions, respectively. We used Poisson statistics in evaluating the upper 
limits on the S*°H*° and f2~f2 _ branching fractions @90CL. The eighth column shows the ratio 
of i/j(2S) to J/ip decays with statistical and systematic uncertainties added in quadrature. The 
last column shows the background subtracted continuum cross-section and the systematic error 
includes 7% due to initial and final state radiation effects and 1% luminosity uncertainty (replacing 
the uncertainty on the number of ip(2S) decays). 
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TABLE II: Breakdown of systematic uncertainty (%) mode by mode. The last column tabulates 
the total systematic uncertainty added in quadrature. 
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In Table HI we show the measured branching fractions of ^(3686) decays to baryon- 
antibaryon modes. The branching fractions for octet baryons are in the range of 0.02 - 
0.035 %. Our measured branching fractions for hyperon modes are ~50% higher than those 
in the PDG |4j which are based on lower statistics (12 and 8 events in the S~ and S° modes, 
respectively). We note that isospin partners, S + and S° and also S~ and 5°, have similar 
branching ratios in agreement with naive expectations. Furthermore, we note that addition 
of strangeness does not greatly change the branching fraction, demonstrating the flavor 
symmetric nature of gluons. The eighth column in Table shows the ratio of these ip(2S) 
results to those from J/ip measurements 0|. Two of the results follow the 12% rule closely, 
but two of them differ by approximately a factor of two. The last column in Table |l] shows 
the background subtracted continuum cross-section at yfs = 3.67 GeV. We quote upper 
limits @90CL for the AA, and production in continuum, owing to mar gina l 

signal. A 20% correction (upward) is included to account for the initial state radiation . 

In conclusion, we have analyzed CLEOIII and CLEO-c ^(3686) data corresponding to 
3.08 xlO 6 ^(3686) decays. We have presented the first observation of the ^(3686) decaying 
to 5°S° and £ + E+ final states, and give improved (high statistics) branching ratios for the 
■^(3686) decaying to pp, AA, and £°£° modes. We also give new upper limits on 

and fi na i s t a tes. 
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